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Abstract 

We consider shape functionals obtained as minima on Sobolev spaces of classical integrals 
having smooth and convex densities, under mixed Dirichlet-Neumann boundary conditions. We 
propose a new approach for the computation of the second order shape derivative of such func¬ 
tionals, yielding a general existence and representation theorem. In particular, we consider the 
p-torsional rigidity functional for p > 2. 
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1 Introduction 


Aim of this paper is to provide a new method for the computation of second order shape derivatives, 
which applies to a broad class of shape functionals associated with classical problems in the Calculus 
of Variations. Let us recall that, if J(-) is a functional depending on a subset R of R n and V is a 
deformation held in G 1 (M ra ;R n ), the first and second order shape derivatives of J at R in direction 
V , if they exist, are given respectively by the limits 


J'(R, V) := lim 
£—^0 


J(R £ ) - J(R) 


J"(R, V) := lim 2 
£—^0 


J(R e ) - J(R) -eJ'(R,V) 


where R £ are the perturbed domains R £ := ^ e (R), with 'f' £ (x) := x + eV(x). 

The literature on shape derivatives is very wide, and in recent years it has seen a rapid flourishing, 
also stimulated by the advances in numerical methods for the determination of optimal shapes. We 
limit ourselves to quote the monographs HZl ESI S2], where a lot of bibliographical references in 
this held can be found. Concerning in particular second order shape derivatives, their computation 
is usually irksome, but often deserves some efforts because the study of their sign allows to detect 
whether a critical shape, namely a domain with vanishing hrst order shape derivative, is actually an 
optimal one. As contributions in this direction, let us mention without any attempt of completeness 
the papers mmmmmmmm- In particular, in [36] Novruzi and Pierre proved a quite 
general and helpful structure result about second order shape derivatives. However, this requires 
the a priori assumption that they exist, and demands that the boundary of the domain R where 
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the derivatives are computed is sufficiently smooth (for the structure of shape derivatives around 
irregular domains, see [29]). Moreover, the linear and bilinear forms which appear in this structure 
result have to be identified each time for the particular functional under study, and this turns out 
to be a delicate task (see the discussion in [25] Section 5.9.4]). 

In this work we propose a new approach to second order shape derivatives, which applies under 
rather mild regularity assumptions to shape functionals of the form 

J(12) := — inf Ij [/(Vu) + g{u)\ dx : u 6 H(fi) j> . (1-1) 

Here / and g are assumed to be convex and smooth functions satisfying suitable growth conditions, 
while H(Vt) indicates the space of functions in H l {£l) which satisfy the Dirichlet condition u = 0 
on a fixed, nonempty, measurable portion T d C <9fl; the Neumann part of the boundary, i.e. the 
complement dQ \Td, will be denoted by IV. 

As a natural continuation of our previous paper [8], where we introduced a new approach to first 
order shape derivatives for functionals of the type m, here we tackle second order shape derivatives 
by the same method. The main results obtained in [8] are briefly recalled in Section [2] see eq. (12.8ft - 
(12.9ft . The basic ingredients we employ are the analysis of the T-convergence of the differential 
quotients appearing in the definition of J"(Q,V), and the duality principle 

J(fi) = J*(9) :=mi{J[f*(a) + g*(dwa)]dx : a € X(Q; R n )j , (1.2) 

where X(f2;R”) is a suitable space of vector fields. 

Our main results, which are intimately related to each other, are: 

- an existence result for J"(Q, V) which is a quadratic form in V and can be represented as follows: 

j"(n,V)= [ (C(u, V) ■ n) dl-L n ~ l + q(u, V ); (1.3) 

Jan 

here C(u , V) ■ n is the normal trace of a suitable vector field (see (13.4ft 1 depending on the solution u 
to J(f4) and quadratically on the deformation V, whereas q(u, V ) is a nonlocal term which involves 
a further vector field B(u, V) and a quadratic form Q(u, •) depending on the second derivatives of 
/ and g (see Theorem I3.4[i : we stress that, as detailed in Remark 13.61 below, formula (11.31) fits the 
general representation result for second order shape derivatives given in [36 , Corollary 2.4]; 

- a regularity result of type for the solution to J(H) (see Proposition 13.2ft : 

- a new necessary optimality condition, which involves the distributional divergence of the above 
mentioned vector field B(u,V) (see Proposition 13.3ft . 

Both the regularity of the solution and the optimality condition on one hand seem to have an 
independent interest, and on the other hand are strictly related to the second order shape derivative. 
Namely, if the regularity of It stated in Proposition 13.21 extends up to the boundary, the field 
C(u, V ) turns out to admit a normal trace on the boundary as soon as the latter is piecewise C 1 ; 
moreover, in order to arrive at the expression (11.3ft for the optimality condition given in 

Proposition 13.31 is exploited as a crucial tool. 

The paper is organized as follows. After providing some notation and preliminary background in 
Section[2l we state in Section[3]our main results lProposition l3.2l Proposition [3T3] and Theorem l3.4l) . 
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followed by some comments and examples. In fact, we show that we can recover straightforward 
the second order shape derivative of the torsional rigidity, as given e.g. in [ 25] Section 5.9.6], and we 
are able to extend the formula to the case of mixed Dirichlet-Neumann conditions. Moreover, we 
consider the case of p-torsional rigidity for p > 2. In this case, to the best of our knowledge, only 
the first order shape derivative is available in the literature (see mm)- We are able to compute 
the second order shape derivative under an additional assumption, which is quite technical and is 
related to a deep regularity issue for the p-Laplace operator; checking its validity seems to be a 
delicate problem of independent interest which goes beyond the scope of this paper. 

The next three sections are devoted to the proofs: the existence of the second order shape deriva¬ 
tive (Theorem 13.41 (i)) is achieved in Section EJ the regularity result of Proposition 13.21 and the 
optimality conditions of Proposition 13.31 are proved in Section [5j eventually, in Section [6] we prove 
the representation formula for the second order shape derivative (Theorem [33 (ii)). 

In Section [7] we present some variants of our results, including in particular the case of the p-torsional 
rigidity for p € [2, +oo), and then we conclude the paper by addressing some possible perspectives. 
Section [8] is an appendix where some auxiliary lemmas and technical facts are collected. 

Acknowledgments. We thank B. Sciunzi for pointing out some useful bibliographical references. 
This work was supported by the University of Toulon and GNAMPA (INDAM). We gratefully 
acknowledge both these institutions. 

2 Preliminaries 

Standing assumptions. Unless otherwise stated, we work under the following hypotheses, which will 
be referred to as standing assumptions: 

- D is an open bounded connected set, with a piecewise C 1 boundary and unit outer normal n; 

/ : R n —>■ R and g : R —>• R are of class C 2 and are strongly convex, namely there exist positive 
constants m and k such that 

(V 2 / — ml) is positively semidefinite, g" — k > 0; 

- / and g satisfy the growth conditions 

Ci(\z\ 2 - 1) < f(z) < C 2 {\z\ 2 + l) R n 

C 3 {\v\ 2 - 1) < g{v) < C 4 {\v\i + 1) WeR, 

where C 7 ; are positive constants, while the exponent q is assumed to satisfy q = 2* := if 
n > 2 and q € (1, +oo) if n < 2; 

- ^(0) = 0 (which is not restrictive up to a translation); 

- the solution u to J(Q) is Lipschitz (we point out that, by strict convexity, u is unique, and 
we refer to HH EH ED E2] for some Lipschitz regularity results); 

- the deformation held V belongs to C 1 (R n ;R n ). 


( 2 . 1 ) 


( 2 . 2 ) 
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Some standard notation. We adopt the convention of repeated indices. Given two vectors a, b in M n 
we use the notation (a, b ) to denote their Euclidean scalar product; moreover we denote by a ® b 
the tensor product of a and b , namely the matrix (a <g> b)ij := afij. Given two matrices A and B 
in W ixn , we write A : B to denote their Euclidean scalar product, namely A : B = AyBy. We 
denote by A~ 1 and A T the inverse and the transpose matrices of A, by A~ T the transposition of 
the inverse of A and by / the identity matrix. Moreover, we denote by a&(A) the k- th invariant of 
the matrix A, in particular 

ai(A)=tr(A), 02 (A) = - (tr(A)7 — A T ) : A , a n (A) = detA. 

For a tensor field A £ C' 1 (M n ; W ixn ), by div A we mean its divergence with respect to lines, namely 
(div A\ := djAij. 

Let ft be as in the standing assumptions. Given a vector field V on <971, we decompose it into a 
tangential and normal component as V = If + V n n, with V n := (V,n). 

We recall that, for ip £ C 1 (0fl) and V £ C 1 (d f2;R n ), the tangential gradient of ip and the tangential 
divergence of V are given respectively by Vr<£> := V</5 — (V<p, n)n and divr V := div U — ( DVn , n) on 
<971, where (p and V are arbitrary C 1 extensions to of p> and V. In particular, if Q is of class C 2 , 
the tangential divergence of n gives the scalar mean curvature, that is Hqq := divp n. The first and 
second order normal derivatives of p> will be denoted by d n (p := {Vip,n} and d 2 n ip := (( V 2 ip)n,n ). 
Given a vector field T £ L 2 (fl;M n ) with distributional divergence in L 2 (Q), we denote by • n its 
normal trace on <90, meant as the unique element in H~ l / 2 (d£i) such that 


{^! -n) ip dK 1 - 1 = 




(2.3) 


where the boundary integral at the l.h.s. is intended as the duality bracket between H 1 / 2 (50) and 
iL 1 / 2 (50). For a detailed account on the theory of weak traces, we refer to [31 H 2 ] (see also [3 
Section 2.1]). 

Existence of solutions, dual formulation and optimality conditions. Under the standing assumptions, 
the inhmum problem J(0) admits a unique solution (c/. [3, Lemma 2.1]), denoted by u, in the space 


H(n) := {u€ H\n) : u = 0onT D }. 


Moreover, by standard duality arguments (c/. [8} Lemma 2.2]), it holds J(0) = J*(Q), where J*(0) 
is the inhmum problem introduced in (11.21) . set over the space 

X(0;M n ) := {a £ L 2 (0;R n ) : diver £ L 9 '(0), a ■ n = 0 on V N } , (2.4) 

being q' := q/(q — 1). Also J*(0) admits a unique solution, denoted by a, which is related to u by 
the following differential equalities: 


a = V/(Vu) a.e. in U 

< 

div <7 = g'(u) a.e. in U 


or equivalently 


Vu = V/*(cr) a.e. in U 

< 

u = (g*)'(div a) a.e. in Q . 


(2.5) 


Note that, by (12.51) and in view of the regularity assumed on /, g and u, we have Vu, a £ L°°(Q] M n ) 
and u, div a £ L°° (f2). 
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First order shape derivative. Let V € C 1 (M n ;M n ) be a deformation field and set for every e 


12 e := \k £ (!2), being \L e (a:) := x + eV(x). 


( 2 . 6 ) 


Under the standing assumptions, the first order shape derivative of J at 12 in direction V, defined 
by 


V ) := lim 
£—^0 


j(n E ) - j(n) 

£ 


(2.7) 


exists and is given by 


j'(n,V)= [ A(u):DVdx= f (A(u)n,V) dU n ~ x , (2.8) 

Jn Jan 

where A(u) denotes the following tensor which turns out to be divergence free on 12 with a normal 
trace A(u)n in L°°(<912): 


A{u) = Vu (g) a — (/(Vu) + g(u)) I = Vu <8> o + (/*(c) + 5 *(div a) — (Vtt, a) — udiv a) I (2.9) 


(with a = V/(Vu)). For the the proof of (12.81) . see Theorem 3.3 and Theorem 3.7 in [SJ; actually, 
even though such results are not stated for the case of mixed Dirichlet-Neumann boundary condi¬ 
tions, it can be easily checked that their proofs continue to work unaltered in such case. Moreover, 
if u is Lipschitz and Vu, a are in SU(12;M ri ), the boundary integral in (12.81) can be rewritten as a 
linear expression of V n , namely it holds J'(Q,V) = h(V n ) with 

h(<p):= [ f\v)vdU n ~ x -\ (/(Vu) + g{u)) dFC 1-1 . 

JTd Jr n 


Reformulation from a variable domain to a fixed domain. Following the same procedure adopted in 
0 Lemma 4.1] in order to recast the first order shape derivative J / (fl, V), we rewrite the variational 
problem J(f2 e ) over the fixed domain 12: by using a standard change of variables, we obtain 

J(n £ ) = - inf / [/ £ (Vu) +g £ (u)\dx (2.10) 

ueH{n) Jn 

where f £ and g £ are given by 

fe(x,z) := f(D^>f T z)/3 £ , g £ (x,v) := g(v)(3 £ (2.11) 

being /3 £ the Jacobian of the map Tg. Note that, in view of (12.61) . f} £ can be written as a polynomial 
in £, with the invariants of DV as coefficients: 

p £ := | det Z2T £ | = 1 + ai (DV)e + a 2 (DV)e 2 + ■ ■ • + a n {DV)e n . (2.12) 

In a similar way, the dual variational problems J*( 12 e ) can be rewritten as 

J*(yt £ )= inf f [/*(<j) + g* (div<r)] dx , (2.13) 

where /* and g* are the Fenchel conjugates of f £ and g £ with respect to the second variable, which 
by direct computation satisfy 

f;(x,z) = r{p- l Dy £ z)P £ , g* E (x,0=g*(fe 1 &vOPe- (2-14) 
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In the subsequent asymptotic analysis as e —» 0, we shall exploit the following developments holding 
for small e, being and f3 e defined respectively by (12.61) and (12.121) : 

D^- t = 1- eDV T + £ 2 {DV T ) 2 + e 3 M £ , with sup £ ||M e || < C , (2.15) 

j3 e = 1 + edivV + e 2 a 2 (DV) + e 3 m e , with sup £ ||m e || < C , (2.16) 

= 1 — edivV + e 2 ((divP) 2 — a 2 {DV)) + £ 3 fn e (V) , with sup £ ||m^|| < C . (2.17) 


3 Main results 


All the results in this Section concern the shape functional J defined in (11.11) and are stated under 
the standing assumptions given in Section [21 Our goal is to provide an existence and representation 
result for the second order shape derivative of J, meant according to the following 

Definition 3.1. Given V € C' 1 (R n ;K n ), J is said to be second order differentiable at in the 
direction V, if the following limit exists: 


J"(n, V) := lim 2 
£—^0 


j(n e ) - j(n) -ej'(n,v) 


(3.1) 


where J 1 is the first order shape derivative of J at Q in direction V (see m> and (BSD for its 
definition and representation formula). 


To prepare our main result, some preliminaries are in order. The next two propositions will be 
necessary in order to prove the main theorem, but also seem to have an autonomous interest: 
as mentioned in the Introduction, the former concerns the regularity of the unique solution u to 
J(f2) (for ir 2,2 -type regularity results for solutions to quasilinear elliptic equations, see for instance 
mmm and references therein); the latter states a new necessary condition for optimality. 

Proposition 3.2. (regularity of the solution) 

Under the standing assumptions, u belongs to W^(fh). 


Proposition 3.3. (necessary optimality condition) 

Under the standing assumptions, the vector field defined by 

B(u, V) := V 2 /(VU)(V 2 U)P - (DV - div V I)a 

satisfies 


div [B(u, V)} = g''(u)(V,\7u) + g'(u) div V in V'(Q). 


(3.2) 

(3.3) 


The vector field B(u,V ) defined in (13.21) will appear in the representation formula for the second 
order shape derivative. Besides B(u,V ), such formula will involve another vector field depending 
on u and V, and a quadratic form on the space i7(fl), which are defined respectively as follows. 
We let C(u, V) be the following vector field, which depends quadratically on V: 


C(u, V) := -(V,Vu)(V 2 f(Vu))(V 2 u)V- (V, Vu)(div VI-DV)Vf(Vu) + 


+ (DV V/(Vu), Vu)V - (DV V, Vu)Vf(Vu) + 


(3.4) 


- (/(Vu) + g(u )) (div VI - DV)V. 
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We emphasize that the distributional Hessian of u appearing in (13.211 and (13.411 belongs to L 1 2 oc (n) 

thanks to Proposition 13.21 

We let Q(u, •) be the quadratic form given by 

Q(u,w) := f 2 [Qf(x,X7w) + Q g (x,w)]dx Vu;ei7(fi), (3.5) 

Jn 

where Qf(x, •) and Q g (x, •) denote the quadratic integrands 

Q f (x,z) := ^ (V 2 f(\/u(x))z,z) Vz G 1", Q g {x,v) := ^ g"(u(x))v 2 VrGl. 

In the sequel for brevity we shall omit to denote the dependence of Qf and of Q g on x and we shall 
simply write Qf(-) and Q g (-) in place of Qf(x, •) and Q g (x, ■). 

We are now in a position to state our main result. 


Theorem 3.4. (existence and representation of the second order shape derivative) 

Under the standing assumptions there hold: 

(i) the second order shape derivative J”(Q,V) defined in (13.111 exists and it is a quadratic form 
in V; 

(ii) ifu is W 2,2 (fiX), the vector fields C(u,V ) and B(u,V), defined respectively in (13.41) and (13.21) . 
admit a normal trace in H~ l / 2 (dUL), and J"(Q,V) is given by 

j”{n,V)= I (C(u, V) ■ n) d% n ~ l + q(u, V ), (3.6) 

Jon 

where q(u,V) is the following non local term: 

q(u,V):=- inf ( Q(u, w - (V, V«)) + 2 f (w - (V, Vw)) B(u, V) ■ nd'H" -1 ) . (3.7) 

w€H( 0 ) 1 Jqq J 


Remark 3.5. Let us show how the representation formula for J"(Q,V ) given by (13.611 can be 
rewritten if one prefers to separate the contributions coming from the Dirichlet and Neumann 
portions of the boundary. First observe that, since u is assumed to be W 2 ' 2 up to the boundary, 
the product (V, Vm) belongs to thus, by considering the traslation v := w — (V , Vu) in (13.71) . 

the nonlocal term q(u,V ) can be recast as 


q(u, V) = — inf 

ueirqn) 

v= — Vnd n u on Fjj 


{ Q(u,v) + 2 [ vB(u,V)-nd'H n ~ 1 } +2 [ (V, Vu) B(u, V) ■ n dli 71-1 
1 h N J JVd 


By combining (13.6|) and (13.81) . we obtain 


(3.8) 


J"(Q, V) = [ ( C(u , V) + 2<V, Vu) B(u, V )) • n dU 71 - 1 + / (C(u, V) • n) dU n ~ l + 

JYb JYm 

inf |q(u, v) + 2 J v B(u,V) ■ nd'H n ~ 1 1; 


(3.9) 


v= — Vnd n u on Fjj 
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in Remark 13.61 below we will see that, under suitable regularity assumptions, the infimum problem 
appearing in (13.911 is a quadratic form in V n . 

Next, in order to further simplify (13.91) . we observe that the boundary integrals appearing therein 
depend only on the normal traces of C(u,V) and B(u,V ) on <917. Therefore, such fields can be 
replaced by any simpler vector field with the same normal trace, which can be chosen in different 
ways on the Dirichlet and Neumann portions of the boundary. On T£>, since Vn is parallel to n on 
<917, C(u, V) + 2(V, Vu) B(u, V) has the same normal trace as 

C D (u,V) := (D,Vu)(V 2 /(Vu))(V 2 7l)D + ((V/(Vn), Vu) -/(Vu))(div V7 - £>V)V; (3.10) 

while on T n, since a ■ n = 0 on such portion of the boundary, C(u, V ) has the same normal trace as 


C N (u, V) := -C V ., Vu)(V 2 /(Vu))(V 2 ?l) V + (V, Vu)DVVf(Vu) + 


+ (DV V/(Vu), Vu)V - (/(VU) + g(u )) (div VI - DV)V. 
We end up with the following reformulation of (13.611 : 


J"{ 17, V ) = / (Cd(«, F) • n) dU n ~ l + f ( C N {u , 1/) • n) dW n_1 + 

Jr jv 

|Q(T, u) + 2 y vB(u, V) ■ ndlT 1 ^ ; 


inf 

v(£H (Q) 

v= — V n dnU on Fjj 


(3.11) 


note that we readily obtain the expression of the second order shape derivative in the pure Dirichlet 
case, i.e. when T^ = (717, by dropping the integrals over T^r in the above equality. 

Remark 3.6. It is interesting to compare our representation formula (13.611 with the structure 
result proved in [SBj, which allows to enlighten the role played by the tangential and the normal 
components of the deformation field V on <917. Notice that, in general, these components are not 
decoupled, i.e. J"(Q,V) / ./"(17, IT) + J"(Q,V n n). In fact, according to the structure theorem 
given in [36l Corollary 2.4, Remark 2.10] (see also [251 Theorem 5.9.2]), under suitable regularity 
conditions J"( 17, V) can be decomposed as J"(17, V) = l\{z) + hiVn), where I2 is a quadratic form 
on C' 1 (917) and l\ is the linear form on (7° (517) associated with the first order shape derivative, 
evaluated in 

z:={V r ,DnVr)-2{Vr,V r V n ); (3.12) 

in this way the scalar 2 encodes the coupling between the tangential and normal components of V. 
In the case of our functional J, by Theorem 3.7 in [8], the linear form l\ is given by 


h(<p) 


f*(a)vdH n - 1 -l lf(Vu) + g{u))ipdU 


'r D 


n —1 




(3.13) 
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and we claim that, as soon as <9fl is C 2 , by Theorem 13.41 the quadratic form I 2 is given by 


h{p>) ■= [ ip 2 [d n u(Dcrn,n) + f*(a)H dn \ dU n x + 

Jr D 

+ f (a, Vri^dnll)) - ip 2 [((V 2 u) W, n) + (f(Vu) + g(u))H m + d n u(DWn, n)] 

Jr N L 

|q(u, v) + 2 J v [ip(Dan, n) — (a, Vr<£>)] . 


<M 


n—1 


+ 


— inf 

v£H^- (fi) 

v= — (pdnu on Fjj 

(3.14) 

This claim can be easily checked when Vp = 0 and <912 = Tpj. Indeed in this case we have z = 0 and 
hence l\{z) = 0; moreover, starting from equality (13.61) . if we replace C(u,V) by the simpler field 
Cd{u,V) given in (13.101) and we use the identity ((divVT — DV)n,n) = V u Hqq (holding thanks 
to the C 2 regularity assumption on <912), we are led to the equality J"(Ct,V ) = hiVn), with I 2 
defined according to (13.141) . The proof of the claim in the general case requires some tedious but 
straightforward computations, that we omit for the sake of conciseness. 

Remark 3.7. We point out that the nonlocal part of the second variation, namely the term q(u. V) 
defined in (13.71) . may also be reformulated in dual form. In analogy with (13.51) . let us introduce the 
quadratic form Q*(cf, •) on the space X(12;M n ') defined in (12.41) : we set 

Q 


*(cr,r]) := / 2[Q f * (x, rj) + Q g * (x, div 77 )] dx , 
Jn 


with Qf*(x,z) := \ (V 2 f*(a(x))z, z) and Q g *(x,v) := \ (g*)"(div a(x))v 2 . (For brevity in the sequel 
we adopt the notation Qf*(-) and Q g *(-) in place of Qj*(x, •) and Q g *(x, •).) 

A standard duality argument (see Lemma 1 8.5 [) yields 

q(u,V) = inf \ Q*(a,rj - B(u,V)) + [ 2 d n uV n (77 • n) d% n_1 ) , 

or equivalently, via a translation, 


q(u, V) = inf 

ri£X(fl\R n )-B(u,V) 


\q*(v,v)+[ 2d n uV n (ri-n)dn n - 1 }+[ 2 d n uV n (B(u, V)-n) dU^ 1 

1 dvn } Jr n 


Example 3.8. An example of shape functional which is covered by Theorem 13.41 is the torsional 
rigidity under mixed Dirichlet-Neumann boundary conditions: 

<7(12) = — inf | j —J-A u^j dx : u £ 77(12) | . 

In this case, assuming 912 of class C 2 and taking deformations normal to the boundary, the second 
order shape derivative given in (13.111) reads 


J' 


\n,v)= v 2 (2\d n u + \d n u\ 2 H m )dH n - 1 + 

z dr D 

[ V 2 [(—2Au + \Vu\ 2 )H dn + 2{V 2 uVu, n)] dU n - 1 + 
^ dr N 

{ [ \Vv\ 2 dx-2 [ v(Wu, VrW) dU n ~ 

L J n J r n 


(3.15) 


inf 

v£H 1 (fi) 

v= — V n d n u on F jj 
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The above equality can be established by applying Theorem 13.41 fas it continues to hold for g(t) = 
—At, cf. the variant discussed in H7.ll) . Namely in this case, as already pointed out in Remark 
13.61 the second order shape derivative is given by fafVn), being I 2 the quadratic form defined in 
([3.141) . The equality between faiV,f) and the r.h.s. of ([3.151) readily follows by using the identity 
An = <9 2 n n + d n uHQn on To, the fact that d n u = 0 on Tat and by taking into account that Q g = 0 

and Qf(z ) = -4-. 

We underline that in the pure Dirichlet case, i.e. when Tjv = 0, our formula (13.151) allows to recover 
the known expression of the second order shape derivative of the torsional rigidity (see for instance 
m Section 5.9.6]). Finally let us point out that, when the deformation field is not normal to the 
boundary, the expression of J"(Q.,V) is slightly more complicated than ([3.151) . Precisely, it is given 
by l\(z) + h(V n )i being l\ and z defined respectively in (13.131) and (13.121) : in other words, we have 
to add to the r.h.s. of formula (13.151) the following terms: 

f ^^-((Vr, Dn Vr)—2(Vr, V r 14)) d'H n - 1 — [ ((Vp, DnV r )-2(V r , V r 14» dU n ~ l 

Jr D 2 j Tn \ z / 


Example 3.9. As a variant of the previous example, we can consider the p-torsion problem under 
Dirichlet boundary conditions: 

J p (D) = — inf | j I-A uj dx : u £ 1T 0 1 ’ P (D)| , (3.16) 

where p > 2 and A is a positive parameter. Note that in this case the integrands / and g do not 
satisfy the standing assumptions, in particular (12.11) . so that we cannot apply directly Theorem 13.41 
In fact a major difficulty appears with the degeneracy at the origin of the Hessian of /. Nevertheless, 
by exploiting a suitable approximation argument, we are able to show that Theorem [23] is still valid, 
provided a suitable equality between weighted Sobolev spaces holds (cf. eq. (17.81) '). The outcoming 
formula reads as follows. Let <9f2 be of class C 2 and let the deformation field V be normal to the 
boundary. Observe that a is parallel to Vlt, and satisfies diver = —A; moreover, the equalities 
((divoT — Da)n,n) = WnHg^ and A u = d^ n u + d n uHgn hold on dQ. Taking these facts into 
account, starting from (13.111) we obtain: 


j^(n,v) 


1 

p 



V*(p\d n u + \d n u\ p H m ) dU n ~ l 


inf / (P(u)Vv, Vv) dx , 

veH qn) Jq 

v= — Vn d n u on df2 


(3.17) 


with 


The complete presentation of this variant is postponed to H7.21 see Theorem O 


Vu 


Vn A 


4 Existence of the second order shape derivative 

This section is devoted to the proof of Theorem 13.41 (i). We introduce for brevity the differential 
quotients 

j(n e ) - j(n) -eJ'(n,v) 


r £ (V) := 2 


(4.1) 
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By exploiting the formulation (12.101) of J(l\) and taking therein as a test function u + ew, we infer 


r £ (V) > ^ ( J [/ e (Vu + eVw) + 5 e (u + ewj] dx - J(Q, V ) - eJ'(fi, L)^ . (4.2) 

We are thus led to introduce the sequence of functionals defined for w £ H(M) and V £ C 1 (M n ; M n ) 
by 

E e (w, V) := |y [/ e (VTH-eV«;)+g E (H + eu;)]dx + J(n)+ eJ'(n, . (4.3) 

Moreover, still for re € i/(fi) and 1/ € C 1 (M n ;M n ), let us define the following functional E(w,V), 
which will be seen to be the variational limit of E e (w, V), and its “dual counterpart” E*(rj, V) (with 
rj £ X(fi;R n ) and V as above): 


E(w,V) := 2 f [f(VTi)+g(u)\a 2 (DV)dx+ 2 \ I [Q f (Vw - DV T Vu) + Q g {w)\ dx + 
Jfi l Jn 

+ J [—((DV — div V" /)V/(Vu), Vtc — DV T 'Vu) dx + div Vg'(u)w\ dxj , 


(4.4) 


E*(rj,V) := 2 / [f*(cr) + g* (diva) — (V/*(cx), a) — u div afiDV) dx+ 


+2 


| j Q f* ((DV — div V I)a + 77 ) + Q g * (div rj — div V div a) + (V f* (a ), DV rj) dx j 


Theorem 13.41 (i) holds in view of the following result: 


(4.5) 


Proposition 4.1. Under the standing assumptions, the second order shape derivative J"(D,V ) 
exists, and is given by 

J"(Q, V) = lim ( — inf EJw,V)) = - inf E(w,V) = inf E*(rj,V). (4.6) 

v ’ e—>0' w&H(n) v ” weH(n) y ’ T,ex(fi;R") w ’ v ' 


Moreover the map V 1 —> J"(D, V ) is a quadratic form in C 1 (M n ;M n ). 


Proposition 14. II will be obtained as a consequence of the two lemmas stated hereafter. 


Lemma 4.2. (upper and lower bounds) 

Under the standing assumptions, the differential quotients defined in m satisfy 


liminf r £ (V) > — inf E(w,V) (4.7) 

e->o ~ w eH ( n) v v ’ 


and 

lim sup r e (fV) < inf E*(rj,V). (4-8) 

e-x) V ex{n-,R”) 


Lemma 4.3. (r-convergence) 

Under the standing assumptions, for every V £ C 1 (M ra ;M n ) the sequence {E e (-,V)} is equicoervice 
and r -converges to E(-,V), with respect to the weak topology of H(Q). 
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Let us show how the proposition follows from these lemmas, and then turn back to their proofs. 

Proof of Proposition \4-l\ We begin bv proving the last equality in (14.61) . Let T : L 2 (14;]R n )xL 9 (n) —>• 
M be the function 

i&(z,v) := / [Qf(z - a) + (z - a, b) + Q g {v) + va + 7 ] dx (4.9) 

Jn 


with 


a:=DV T 'Vu, b := — (DV — div V I)a , a := divVg'fu), 7 := [/(Vn) + g(u)] a 2 (DV). 

For every w G it holds E(w, V) = 2 '&(S7w,w), in particular 

— inf E(w, V) = — 2 inf {^(\/w,w)} . (4-10) 

weH(n) w£H(Q) 

By applying Lemma 18.51 with Y = H(Q), Z = L 2 (fi;M n ) x L q (Q), A : Y —> Z the operator 
Av := (Vr,r), $ : Y —> R the zero function and 'F defined as in (14.9p . we can rewrite (|4.1(JI) as 

— 2 inf {^(Au;) + 4>(rc)} = 2 inf {'F*(ry, r) + A*(rj, r))} . (4.11) 

w&Y (v,r)ez* 

Let us compute the Fenchel conjugates appearing in the r.h.s. of (14.111) . By exploiting Lemma 18.21 
it is easy to see that 

^*(V, r) = [ [( Qf)*(v ~b) + {tj, a) + ( Q g )*(r - a) + 7 ] dx . 

Jn 

Since <f> = 0, its Fenchel conjugate < F* is 0 at 0 and +00 otherwise. As an element of Y*, A*(rj,T ) 
is characterized by its action on the elements of Y: since 

{A*(ti,t),v)y*,y = Av) z *,z = / (( 7 , \7v)+tv) dx , 

J n 

we infer that <&*(A*(rj, r)) = 0 if and only if r = div 7 / (with the additional condition 77 • n = 0 on 
Tjv); note that, as a consequence, such vector fields r] belong to A(fl;M ri ). 

Therefore we may rewrite the r.h.s. of (14.111) as 

2 Jjo [ KQ fY (V ~ b) + (v, a) + (Q g )* (div r)-a)+ 7 } dx) . (4.12) 

77ex(n ; R-) (.Jn J 

Since by Lemmas 18.31 and 18.41 we have ( Qf)* = Qf* and (Q g )* = Q g *, and by the optimality 
conditions (12.51) there holds 

7= [ [f*(a) + g*(div a) - (Vf*(a), a)-u divw]a 2 (DV)dx , 

Jn 

we conclude that the expression in (|4.121) agrees with inf^g wq ; r») E*(i 7 , V ), which achieves the proof 
of the last equality in (14.61) . Combined with Lemma 14.21 such equality implies the existence of the 
second order shape derivative, and the fact that it agrees with any of the two infimum problems at 
the r.h.s. of (14.71) and (14.81) . The second equality in (14.61) is a direct consequence of Lemma FTTol 
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Finally, the fact that the map V H > J"(Q,V) is a quadratic form is readily checked: one may 
start from the equality J"(Q, V) = — inf,,. E(w, V), and look at the expression 1)4.41) of E(w, V ), by 
recalling in particular that Qf is a quadratic form. □ 


Proof of Lemma \4-°2\ We first prove the lower bound in (14.71) . We start from the inequality (14.21) . 
where we choose w as an element of C°°(fl) fl H(Ll). Recalling the definitions (12.111) of f £ and g £ , 
the definition ( 11 . 11 ) of J(fl), and the expression (12.81) of V), we can rewrite the inequality (14.21) 


separating the terms of different orders in e. 
according to (12.151) and (12.161) . and setting for brevity 


Recalling that M e and m £ are defined respectively 


hi :=Vw-DV t Vu, h 2 := (DV T ) 2 Vu - DV J Vw , h £ := (DV 1 ) 2 Vw + M e (Vu + eVw) 


rT\ 2v7tt 


T\2\ 


we get 

r £ (V) > (/ 0 (e) + h{e) + hie)) , (4.13) 

where 

Iq(s) := — 2 [/(Vu + eh\ + e 2 h 2 + £ 3 h £ ) + g(u + ew)](a 2 (DV) + em £ ) dx , 

Jn 

/i(e) :=- f [/(Vu + eh\ + e 2 h 2 + e 3 h £ ) + g(u + ew) — /(Vu) — g(u)\ div V dx , 

e Jn 

h{e) ■= ~ ~o [ [/(Vu + eh\ + e 2 h 2 + e 3 h £ ) + g(u + ew) - /(Vu) - g(u) + e(V/(Vu), DV T 'Vu)} dx . 

£ Jn 

Let us study separately the asymptotic behavior of //e) as e —>• 0. By exploiting the growth 
assumptions on the integrands and the uniform L°° boundedness of m £ and h £ , we get 


lim J 0 (e) = -2 f [/(Vu) + g(u)]a 2 (DV) dx . (4.14) 

By applying Lemma 18.11 (i) to the integral functionals If and I g defined according to (18.11) (note 
that this can be done after possibly modifying g outside an interval, since u + ew remains uniformly 
bounded in L°°), we get 


lim Ji(e) = -2 / [(V/(Vu), Vu; - DV T \7u) + g\u)w] div V dx . (4.15) 

£ ->° Jn 

Let us now consider I 2 (e). We recall that, in view of the optimality conditions in (12.51) . there holds 
a = V/(Vu) and diva = g'(u ), so that 


[(V/(Vu), X7w) + g'(u)w] dx = / ((a, Vw) + w div a) dx = 
i Jn 


wo 


ndU n ~ 1 = 0, (4.16) 


'an 


where the last equality holds since w = 0 on T/r and o-n = 0 on r^r. Thus, I 2 (e) remains unchanged 
if we add the zero term (14.161) (multiplied by 2/e): 


hie) 


-o / [/(Vu + e/ii + e z h 2 + s d h £ ) — /(Vu) 
e Jn 

-o / [g(u + £ w) - g{u) - eg\u)w} dx . 

£ Jn 


e(V/(Vu), h\)} dx + 
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Hence, by applying Lemma l8.il we obtain 


lim/ 2 (e) = -2 / [Q/(M + (V/(Vu), h 2 ) + Q g {w)] dx . (4.17) 

Jn 

By combining (14.131) . (14.141) . (14.151) and (14.171) . by the arbitrariness of re € D H(fl) we infer 

lirninf r £ (V) > — inf E(w,V ). (4-18) 

£ ^° weC°°(n)nH(n) 

In order to conclude the proof of the inequality (14.71) . it is enough to observe that the inhmum at 
the r.h.s. of (I4.18P coincides with the inhmum at the r.h.s. of (14.71) . Indeed, the integral functional 
E(-,V) above is continuous on H(Q) (by the growth conditions (12.21) and the Sobolev embedding 
of W 1 ,p (£l) into L p ( H)), whereas C' 0 O (H) n H(£i) is dense in iL(fi). 

Let us now prove the upper bound in (14.81) . Let g be an arbitrary element of the space C°°(H; M n ) n 
X(fl;R”). In view of (11.21) . the differential quotients r e (V) introduced in (14.11) can be rewritten as 

r e (v) = J, [J*m - r(n)-ej'(n,v)] . 

Then, by exploiting the expression of J*(fi e ) in (12.131) . we infer that 

r £ (V) < ^ (^J [f*(a + £g) + g*(diva+ £divg)] dx - J*(Q) - £j'(Q,V)) . (4.19) 

Recalling the definitions (12.141) of /* and g *, the definition of J*(H) in (11.21) . and the second 
expression for J'(Q] V ) in (12.81) . we can rewrite the inequality (14.191) separating the terms of different 
orders in e. Recalling the definition (12.171) of fn £ , and setting for brevity 

zi :=(DV — di vVI)a + g , 

z 2 =(DV - div V I)g + ((div V ) 2 - div V DV — a 2 (DV))a , 

z £ :=fn £ (W + £rj) + ((div V ) 2 — a 2 (DV))(DVa + g) — divR DVg , 

r\ := div g — div V div a , 

t 2 :=((div V ) 2 — a 2 (DV )) div a — div V divg , 

r £ :=((div V ) 2 — a 2 (DV)) divg + m £ (div u + ediv g ), 


we obtain 

re{V) < (/q( e) + It{e) + 12(e)) , (4.20) 

where 

Iq(£) :=2 f [f*(a + £Zi + £ 2 z 2 + £ 3 z £ ) + g*(a + £T± + £ 2 t 2 + £ 3 T £ )\(a 2 (DV) + £m e ) dx , 

Jn 

I* (e) :=- [ [/*(u + ezi + £ 2 z 2 + £ 3 z e ) + g* (a + £T\ + e 2 r 2 + £ 3 t £ ) - f*(a) - g*(diva)] divR dx , 

£ Jn 

I 2 (e) : =4 f [f*(v + £Z 1 + + £ 3 z £ ) + g*(a + ETi + £ 2 T 2 + £ 3 T £ ) - f*(o) - g*(diva) + 

e 2 Jn 

— e(V/*(cf), ( DV — div V)a) + ediv V div a(g*)' (diva)] dx . 
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Let us study separately the asymptotic behavior of I*(s) as e —> 0. By exploiting the growth 
properties of the integrands (note that the growth assumptions on / and g made in ( 12.211 imply 
similar conditions on /* and g*) and the uniform L°° boundedness of fn £ , we get 

lim/g(e) = 2 [ [f*(a) + g*(diva)]a 2 (DV) dx . (4-21) 

£ ->° Jn 

By applying Lemma [8711 (i) to the integral functionals If* and I g * defined according to () 8 . 1 f) (note 
that this can be done after possibly modifying g* outside an interval, since diva + edivg remains 
uniformly bounded in L°°), we get 

limine) = 2 f [(V/*(a), zf) + (g*/(div a)ri] div V dx . (4.22) 

£ ->- 0 J n 

Let us now consider /|(e). We recall that, in view of the optimality conditions (|2.5j) . we have 


/ [(V/*(<x),r/) + ( 3 */(diver) div 77 ] dz = 0. (4.23) 

Jn 

Thus, 7|(e) remains unchanged by adding the zero term (14.2311 (multiplied by —2/e): 

Ifa) =~2 [ [f*{° + ez i + e 2 z 2 + £ 3 Zc) ~ f*(&) -e(Vf*(a),zi)+ 

£ Jn 

+ g*(a + £T\ + e 2 T 2 + e 3 r e ) — 5 f*(div a) — e( 5 f*) , (diva)ri] dx . 

Hence, by applying Lemma 18.11 we infer 

lim/^e) = 2 f [Q f *(z 1 ) + (' Vf*(a),z 2 ) + Q g *(n) + (g*)'(div a) r 2 ]dx. (4.24) 

Jn 

By combining (14.2011 with (|4.21l) . (I4.22P and (|4.24|) and recalling the definition (14.511 of E*, in view 
of the arbitrariness of g € C' 00 ^;!!”) n 7T(f2;M n ), we infer 

limsupr e (V r ) < _ inf E*(g,V). (4-25) 

e-i-0 r?eC 00 (r2;K n )nX(r2;R 11 ) 

Finally we observe that, for every g £ X(H;M. n ), there exists a sequence {%} C C°°(n;R n ) n 
X(ri;]R n ), such that r/h —> g in L 2 (14;M n ) and di vgh — > divg in L q '(U). Since E*(-, V) is continuous 
with respect to such convergence (thanks to the fact that the growth assumptions on / and g made 
in (12.211 imply similar conditions on f* and g*), the infimum at the r.h.s. of (I4.25H turns out to agree 
with the infimum at the r.h.s. of (14.811 . and the proof is achieved. 

□ 


Proof of Lemma f.3. For brevity, throughout the proof we fix the vector field V and we omit it 
when writing the functionals E e and E. We begin by showing the equicoercivity of E £ . In view of 
the convexity assumptions on / and g we infer that the following bounds hold true, where m and 
k are the positive constants appearing in ( 12 . 11 ) : 


f(z 0 + z) > f(z 0 ) + (Vf(z 0 ),z) + m\z \ 2 , g(v 0 + v) > g(v 0 ) + g'(v 0 )v + k\v \ 2 . 
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By exploiting such bounds for zo = Vu, z = —Vu + D\ L e T (Vu + Vw), v$ = u and v = sw, we infer 
that for every e 


E e {w) > ci ||Vu>||^ 2 + c 2 |Ml| 2 -c 3 \\Vw\\ L 2 -c 4 |M| L 2 - c 5 , 

for some positive constants c*. Therefore the sequence E £ is equicoercive in H(H). 

Let us now show the T-convergence statement. By definition of T-convergence, we have to prove 
the so-called T-liminf and T-limsup inequalities: 



<| lirri inf E £ (w £ ) 
lim sup ^(u^) 


H(n) 

V £ V 


H(fl) 

V £ V 


■ > E(w) 

(4.26) 

< E(w). 

(4.27) 


Let us prove (14.261) . For every w £ H(Q) and for every sequence w £ which converges weakly to w 
in H(Q), we have to prove that 

E(w) < liminf EJw £ ) . (4.28) 

£ 

For a given w € H(Q), let us first consider a regular sequence w £ £ C°°(£l)nH(Q) weakly converging 
to w. 

As already done in the proof of Lemma r4.2l . by exploiting the C 1 regularity of the integral functionals 
If and I g and the L 2 -Mosco convergence of the corresponding sequences A £ j and A £;g dehned 
according to (18.21) . we infer 

E(w) < liminf E £ (w £ ). (4.29) 

£ 

Let us now consider a generic sequence w £ £ H(Q) weakly converging to w. By density, for every 
fixed e, there exists w £) k £ C°°(n) D H(Q) such that w £: k —> w £ strongly in H(Q) as k —> +oo. With 
a diagonal argument, we may find a subsequence such that 

lim inf E £ (w £ ) = lim inf lim E £ {w £ *.) = lim EJw e u) , (4.30) 

e £ k ’ £ ’ 

where the first equality follows by the strong continuity of E e . We remark that the subsequence 
w £y k e is regular and weakly converging to w, hence it satisfies (|4.29l) . Such property, combined with 
the equality (14.301) . concludes the proof of (14.281) . 

Let us now prove (|4.27j) . For every fixed w £ H(H) we have to find a recovery sequence, namely 
a sequence w £ which converges weakly to w in H(Q) and satisfies limsup^g E e {w e ) < E(w). If 
w is an element of the space C°°(n) D iL(fl), we are done simply by taking the constant sequence 
w £ = w. Indeed, by following the same procedure adopted in the proof of Lemma 14.21 we infer 
that lim e _j.o E e (w) = E(w). If re is a generic element of H(Q), we approximate it by a sequence 
~ g C°°(Q) nH(fl): by the lower semicontinuity of the l.h.s. of (|4.27l) (usually called T-limsup E e ) 
and the continuity of E, we obtain 

(r — lim sup E e )(w) < liminf(r — limsupF/ £ )(?Ufc) < liminf lim sup E e (wk) = E(w ). 

k k e 


□ 
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5 Regularity of the solution and necessary optimality condition 


In this section we prove Propositions 13.21 and 13.31 Both of them are based on the crucial result 
given below, which in turn exploits the T-convergence statement given in the previous section (c/. 
Lemma 14.31) . 


Proposition 5.1. Under the standing assumptions, if V has compact support contained into fl, 

then the function 9y := (V, VU) solves the minimization problem inf E(-,V), where E(-,V) is 

w€H(Q) 


the functional defined by (HI- 


Proof. We first prove the following claim: for every fixed e > 0, the function 

uix + sVfx)) — u(x) 
w £ {x) :=- : - 


(5.1) 


solves the minimization problem inf E e (-,V), where E e (-,V) is the functional defined by (14.31) . 

w&H(Q) 

Let u e G H{Ul) be optimal for problem J{Ul £ ) in the formulation (12.101) set on the fixed domain fh 
Via change of variables we obtain 


J(fl £ ) = - f [/ £ (Vu £ ) +g £ (u £ )]dx = - [ [/(V(it e o 1 )) + g{u £ o 1 )] dx . 

We infer that the function ri £ o'I '“ 1 € H(Q e ) is optimal for problem J(O e ) in its original formulation 
as an infimum over Ef(Q e ). On the other hand, since V has compact support contained into 14, we 
know that J(I4 £ ) = J(I4), and hence (since u is the unique solution to J(I4)) we have u £ o T ” 1 = u, 
namely u e (x) = u of e (x) = u{x + eV{x)). It follows that the function w £ defined in (15.11) minimizes 
E e (-,V) over and our claim on the optimality of w £ is proved. 

In view of Lemma 14.31 the minimizing sequence {tr £ } is equibounded in Ed(Ul) and, up to sub¬ 
sequences, it weakly converges to a minimizer of the T-limit, namely to a minimizer of E(-,V). 
Clearly, the weak limit of w £ necessarily agrees with the pointwise limit (V, Vn), which allows to 
conclude that (V, Vu) minimizes E(-,V) over H( fl). □ 


Proof of Proposition \ In view of Proposition 15.11 we can assert that, for every vector field 

V compactly supported into P, the function (V. Vu) belongs to By the arbitrariness of 

V € C ( J (14; M”), we infer that u is an element of □ 

Proof of Proposition [PI We consider the minimization of the functional E(-, V ) in (14.41) over 
written as done in the proof of Proposition 14.11 see eq. (14.101) . By duality, we have the equality 
(14. 11|) . Then, by Lemma 1^31 a function w G H(Q) and a pair ( 77 , r) G L 2 ( fl; W l ) xL q (14) are optimal 
respectively for the two problems in (14.111) if and only if (7/, r) G d^CVw, w) and (div 77 , r) G d&(w). 
Then, by computing explicitly the expressions of <9$ and cM, we get r = div 77 , with the additional 
condition 77 • n = 0 on Tjv, and 


77 = V 2 /(Vm)(Vtc - DV t Vu ) - (DV - di vVI)a, 

< 

div 77 = g"(u)w + g'(u) div V . 
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Hence, for V £ Cq(H;M”), the distributional equality (13.31) follows by replacing w by (V, Vu) in 
the above optimality conditions (since we know by Proposition 15.11 that 9y = (V, Vm) minimizes 
E(-,V) over H(Q)). For V £ C' 1 (R n ;M n ), the equality (13.31) continues to hold by the following 
simple argument: in order to check (13.31) for a given test function <p £ X>(P), it is enough to use its 
validity when V is replaced by the product gV £ C'q(P;R ? 1 ), being g a cut off function compactly 
supported into H which equals 1 on the support of p. 

□ 


6 Representation formula for the second order shape derivative 

In this section we give the 

Proof of Theorem \3.Jj\ (ii). We start from the following equality that we have already proved in 
Proposition Id. II 

J"(n,V) = - inf E(w, V). (6.1) 

Throughout the proof we set for brevity 9y := (V, Vu). By inserting into the definition (14.41) of 
E(w,V) the identities 


Q f (Vw - DV t Vu ) = Q f (V(w - 9 V )) + (V 2 /(Vu)(V%)P, V(w - 9 V )> + Q f ((V 2 u)V) 

Qg(w ) = Qg(w- 9y) + g"(u) 9y (W - 9y) + Q g (9y) , 
we obtain 

E(w,V) = [ 2 [Q f (V(w-0y))+Q g (w-e v )] dx + 


+ / 2[(V 2 /(Vu)(V 2 u)P- (DV - dwVI)a,V(w - 9 V )) + (g"{u) 9y + divVg'(u))(w - 9 V )] dx + 


+ 


[ [(/(Vu) + g(u)) (div VI - DV) : DV T - 2{(DV - div VI)W, (V 2 u)P) + 2 div Vg'(u)9y\ dx + 
Jn 


+ / [2Q f ({X7 2 u)V) + 2Qg(9y)\ dx. 

Jn 

Recalling the definition (13.51) of the quadratic form Q(u, ■) and the definition (13.21) of the vector field 
B(u, V), by using also Proposition 13.31 we can rewrite the above equality as E(w,V ) = (I) + (II), 
with 

(I ) = Q(u,w — 9y) dx + [ 2[(B(u, V), V(w — 9y)) + divB(u, V)(w — 9y)\ dx 

Jn 

and 

(II) = [ [(/(Vu) + g(u)) (div VI - DV) : DV T - 2{(DV - div VI)a, (V 2 u)P) + 2 div Vg'(u)9 v \ dx + 
Jn 

+ [ [2Q f ((V 2 u)V) + 2 Q g (9y)\ dx . 

Jn 
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Note that, in view of the assumption u £ W 2 ’ 2 (Q), the vector field B(u, V) is in L 2 and has bounded 
divergence, thus it admits a normal trace in H~ l ^ 2 (dV), and we are allowed to apply the integration 
by parts formula (12.31) . Thus we get 

(I) = Q(u, w -e v )dx+ [ 2 (w-Oy)B(u,V) ■nd'H n ~ 1 , ( 6 . 2 ) 

Jan 

We now want to rewrite (II) as a boundary integral. To this purpose we observe that, by (13.21) and 
(J3.3D . there hold 


-((DV - div VI)a — V 2 /(V 2 u)V, (V 2 u)U) = ( B(u,V ), (V 2 u) V) , (6.3) 

div Vg'(u ) By + g"(u)(By) 2 = div B(u, V) By . (6-4) 

Moreover, again thanks to the assumption u £ IT 2 , 2 (II), the vector field 

Xi := (/(Vm) + g(u))(divVI - DV)V (6.5) 


is bounded and has L 2 divergence. Thus we are allowed to apply the integration by parts formula 
(J2.3D and we get 


f (f(Vu)+g(u))(divVI - DV) : DV T dx = [ (/(Vu) + g(u))((divVI - DV)V,n) dH n ~ l + 
isi Jdn 

- [ ((di vVI - DV)V,(y 2 u)a + g'(u)Xu)dx. ( 6 . 6 ) 

Jn 

Using (16.31) . (16.41) and (16.61) . we infer that 


(II) =[ (X! ■n)<m n - 1 + f [(B(u,V),(V 2 u)V) + div B(u,V)9 v ]dx + 

Jan Jn 

- f {(X 2 u)DVa — DV T (y 2 u)a — g' (u)DV T 'Vu, V) dx . 

Jn 


Exploiting the equality (V 2 u)V = V6y — DV T Xu , we get 


(II) = (Xi •n) < m n " 1 + / [(B(u,V),XB v ) + div B(u,V) By] dx + 
Jan Jn 


in L 


(V 2 /(Vu)(V 2 u)U, DV t Vu) + div V(a, DV T Xu) + (( V 2 u)DVa , V) 


dx + 


+ 


in L 


(DVa, DV t Xu) + (DV 1 (V 2 u)a, V ) + gf(u)(DV T Vu, V) 


r T/-n2—\— 


! Orr\ / 711 rT V7TT 


dx. 


Finally we remark that (B(u, V), XBy) + di v B(u, V ) By = div ( 9yB(u , V)), and 


- (V 2 /(Vu)(V 2 u)V, DV t Xu) - div V(a, DV T Xu) - (( X 2 u)DVa , U) + (DVa, DV T Xu)+ 
+ (DV T (X 2 u)a, V) + </(u)(.DU r Vu, U) = div ( - (DVa, Vu)V + (DV V, Vu)a ) . 
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Hence 


(II) = [ (X 1 -n)dH n ~ 1 + [ div (d v B(u,V) - {DVa,Vu)V + (DVV,Vu)a) dx ; (6.7) 

Jan Jn ' ' 

moreover, thanks to the assumption u £ W 2,2 (Q), the tensor field 

X 2 := 0 V B(u, V) - {DVa, Vu)V + (DV V, Vu)a , (6.8) 

is in L 2 and has L 2 divergence. Then, by adding (|6.2p and (16.7|) and by applying the integration by 
parts formula (12.31) . we obtain 


E(w,V) = Q(u,w — 6y) dx + f 2(w — 0y) B(u, V) ■ n dT-L n x + 

Jan 


+ [ (X, • n) dU n - 1 + f (X 2 ■ n) dU n ~ x 

Jan Jan 


(6.9) 


Now, the fact that the field C(u, V) defined in (13.4|) admits a normal trace in i7 _ 1 / 2 («9H) follows by 
observing that C(u, V) = —X\ — X 2 , with X\ and X 2 defined respectively in (16.51) and (| 6 . 8 I) . Finally 
the expression (13.61) of the second order shape derivative follows by combining (16.91) with (16.11) , and 
using the equality C(u, V) = —X\ — X 2 . 

□ 


7 Variants and perspectives 

7.1 The case g linear 

We claim that Theorem 13.41 Proposition 13.21 and Proposition 13.31 remain true in the case when 
g(v) = — Xv, for some A £ R. Note just that in the definition (13.51) of the quadratic form Q(u, ■), 
the term Q g vanishes. 

The proof of the lower bound inequality (14.71) works unaltered, whereas the proof of the upper 
bound inequality (14.81) must be modified as follows. Note firstly that, since the Fenchel conjugate 
of g is the Dirac delta function g*{r) = 5_\, the dual energy E*(g, V) introduced in (14.5p is finite 
only if div ?7 = — AdivP, and in this case it reads: 


E- 


(v,V)= 2 [ if*(a) - (V/*(<f),<f> + Xu] a 2 (DV)dx + 

J n 

+2 / [Q f *((DV-divV I)a + r]) +{Xf*(a), DVg)] dx . 
Jn 


(7.1) 


Next, we write the following inequality (which replaces (14.191) 1: 

2 


r £ (V) < f 2 ^JJf:(a £ ) + g* £ (dWa £ )}dx-J*(n,V)-£j'(n,V) 

where a e are the perturbations of a defined by 

n 

+^2 £kr]k ’ 


(j, := a 


(7.2) 


(7.3) 


k=1 
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with for k £ { 1 , 2 ,n} 

r/k £ L°°(£l; R n ), div rj^ = —Xak(DV) in Q , g k ■ n = 0 on Tn . (7-4) 

Notice that, for such a choice of a £ , the r.h.s. of the inequality (17.2ft is finite. Namely, exploiting the 
expansion ( 12 . 121 ) and the fact that by (12.5ft it holds diva = —A, we obtain diva £ = — Xf5 e , so that 

g* (div cr e ) = fi £ g* (/3' 1 div a £ ) = (5 e g*(-X) = 0 . (7.5) 

Notice also that the set of fields satisfying (17.41) is non empty: indeed, every a k (DV) can be expressed 
as the divergence of a suitable regular vector field u k (see for instance [33l Lemma 4.6.4]); thus the 
conditions in (17.41) are verified by taking g k := — Xu^ + Vu;*,, being w k the solution of 


A w k = 0 in , d n Wk = h k on , 


with hk = Xuk ■ n on Tjv , f 9Q hk = 0. By combining (17.21) and (17.51) . we may write 

r £ (V)< (/ 0 *(e)+/i*(e) + i!(e)), 


(7.6) 


with 


r 0 (s) :=2 [ f*(p- 1 D^ £ a £ )(a 2 (DV) + em £ )dx, 

Jn 

im ■■=- [[np^DVe^-nn&vvdx, 
e Jn 

I 2 (e) :=4 f [riPe'D^eVe) - /*(<f) - e(Vf*(a), (DV - di vVI)a) - eAdivFu] dx . 

£ Jn 

Now we argue in a similar way as in the strictly convex case: we apply Lemma 18.11 to the integral 
functional If* and, in order to deal with I 2 (e), we exploit the equality 

/ Adi vVudx = — / udivr]idx= / (V/*(a), r/i) dx ; 

Jn Jn Jo, 


thus we obtain 


= 2 

[ f*(cr)a 2 (DV)dx , 

£->0 

In 

lim/r(e) = 2 

[ (V/*(cf), (DV — div V I)a + r/i) div V dx 

£ — ^0 

In 

hni I 2 (e) = 2 
£->■0 

f Qf*((DV — div V I)a + 771 )+ 

ln 


+ <V/*(er), ((div V) 2 - div V DV - a 2 (DV))a + (DV - div V I) m + r/ 2 ) dx . 


In view of (17.61) . by adding up the three terms above and exploiting the equality f n (V f* (a), g 2 ) dx = 
Xua 2 (DV) dx, we infer that limsup^g r £ (V) < E*(r] i, V). 

Now, since rji is an arbitrary vector field in L°°(fl;M n ) with divryi = —A and r/i • n = 0 on T^v, in 
order to conclude the proof of the upper bound inequality (14.8ft it is enough to show that every field 
g £ A(f2;M n ) satisfying the two conditions div?? = —A and rj ■ n = 0 on can be approximated 
strongly in L 2 by a sequence of fields rj h belonging to L°°(f2;M n ) and satisfying the same two 
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conditions. To that aim we begin by noticing that, since by assumption 30 is piecewise C 1 , there 
exists a sequence of smooth fields fj h , with rj 1 ■ n = 0 on Tjv, such that fj h —>• rj in L 2 (0;M n ) and 
di vfj h —> — A in L 2 (0). Then, we let (j) h be the unique solution to the Dirichlet-Neumann problem 

— A(j> h = p h := A + div fj h in 0 , <p h = 0 on T^ , d n (j) h = 0 on T^ ■ 

Finally we define rj 1 := fj h + V(f h , and we claim that such sequence rj h has the required properties. 
Namely, since p h is a sequence of smooth functions converging to 0 in L 2 (0), the sequence q h lies in 
L°°(0;M n ), and converges to rj in L 2 (0;R ra ); moreover it holds divr/' = —A and r] h ■ n = 0 on Tat. 
Eventually, once now (14.81) is established, the remaining parts of the proofs given in Sections [3] [5] 
and [6] can be repeated as done for g strictly convex. 

7.2 The p-torsion problem 

Here we prove formula (13.171) for the second order shape derivative of the p-torsional rigidity func¬ 
tional J P {Q) introduced in (13.161) . for p > 2. We recall that the first order derivative J' p {Ll, V ) exists 
and is given by 


4(n,v)= f ^^v n( m n -\ 

p Jan P 

where u is the unique solution in Wq' p (fi) to —A p u = A (see for instance [ 8 , HU). It is well-known 
that u is of class C' 1 ,Q (fl) (see [18, 33])) and it is also in C 2 (fi \ S ), where S is the critical set 
S := {x £ Ll : Vu = 0}. Let us also recall that S has vanishing Lebesgue measure (see [30]4. and 
it is compactly contained into fl (thanks to Hopf boundary lemma). 

Let us introduce the weight function p := |Vu| p-2 , which is continuous in 0 and strictly positive 
outside S. Moreover, p _1 is in L 1 (0) (see [30]'). 

We denote by kFp’ 2 (fl) the Hilbert space consisting of functions v € H 7 ^ 1 (Q) such that 

ll v llwp' 2 (Q) := ll w lli 2 (n) + H Vv llL2(n) < +°°i 

and by Hl’ 2 (Ll) the completion of C , 1 (fl) with respect to the above norm. 

Let Q(u, •) denote the quadratic functional defined on W p ' (H) by 

Q(u, v) := J (.P(u)Vv , Vu) dx , with P(u) := |Vu| p ~ 2 ^1 + (p - 2 )|^=| ® ■ ( 7 - 7 ) 

Notice that, since the norm of the matrix P(u) is controlled by p , the functional Q(u, ■) is continuous 
on the Hilbert space LFp’ 2 (H). By adapting the approach developed in Section we obtain the 
following second order differentiability result. 

Theorem 7.1. Let dLl be of class C 2 . Let p > 2 and assume that the weight p := |Vu| p-2 is such 
that 

H 1 p 2 (n) = w}’ 2 (n). (7.8) 

Then, the functional J P (Q) is twice differentiable at Ll in any direction V, and for V normal to the 
boundary it holds 

J"(n,V) = --[ V%(p\d n u + \d n u\ p H m ) dli n - 1 - inf Q(u,v). (7.9) 

P Jan n) 

v= — Vndn~ on dQ 
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1 2 

Remark 7.2. Notice that, since p is strictly positive in a neighborhood of dQ, the space W p ’ (fl) 
is embedded into \ S ), and hence the trace operator is well defined from W p ’ (Q) into 

i7 1 / 2 (i90). In particular, the infimum in (17.91) is well-defined since, by the C 2 regularity of u in 
Q\S, the function (V, Vu) belongs to H\ oc (Q\S). Moreover, the infimum can be equivalently taken 
in I7 1 (n) (as written in Example 13.91) . 

Remark 7.3. Let us add a few comments on assumption (17.81) and related bibliographical references. 
As general texts on weighted Sobolev spaces, we refer to [28||33j. Moreover, we address to the paper 
|13| for a counterexample showing that the L 1 -summability property of p -1 is not strong enough to 
ensure the validity of (17.81) . A sufficient condition would be p £ M 2 , where M 2 is the Muckenhoupt 
class of functions satisfying sup B Jb p) Ib P _1 ) < +°°> where B varies among balls in R n 
(see for instance ESI), but such kind of regularity result seems hard to prove. On the other hand, 
equality (17.81) holds true if the set S has vanishing capacity in H P , 2 (D,), namely if there exists a 
sequence of smooth functions a e : 0 —> [0,1] which are equal to 1 in a e-neighborhood of S and 
satisfy lim e f n \\7a £ \ 2 pdx = 0. Indeed in this case, given w £ IEp’ 2 (0) with w = —9y on dii, it 
is possible to construct an approximating sequence {w e } £ H l p ' 2 {D) with w e = — 6 y on dQ: it is 
enough to take w £ := (1 — a £ )w (which belong to L/Rfi), and hence to H p ' (0)). In particular, our 
assumption (17.81) turns out to be satisfied whenever S is a singleton. Concerning the geometry of 
S recall that, if II is convex, the function u is power concave [39], and hence S agrees with the set 
where u assumes its maximum. If in addition II is strictly convex, it is likely true that S is reduced 
to a singleton (c/. [5] 122]), or at least that it is of vanishing capacity, so that (17.81) holds true. 

Proof of Theorem \7.1\ Since shape derivatives only depend on the behavior of the deformation field 
on the boundary, with no loss of generality we may choose V vanishing in a neighborhood of S 
(recall that S CC II). Then, in order to show the equality (17.91) . we proceed along the same scheme 
adopted in Section UJ namely we show that the lower and upper limits as e —> 0 of the differential 
quotients r £ (V) defined in (14.11) are bounded respectively from below and from above by the same 
quantity, which is precisely the r.h.s. of formula (17.91) . For convenience, we divide the remaining of 
the proof in three steps. 

Step 1 (lower bound): Thanks to the assumption p > 2, the following lower bound can be achieved 
by arguing exactly as done in the proof the inequality (14.71) in Lemma 14.21 

lirninf rJV) > moo := — inf E(w, V ), 
e—>-0 v ’ - w£Cg°(n) 

where E(w, V) is defined according to (14.41) by taking therein f(z) = \z\ p /p, and g(v) = —An. 

We claim that the field B(u,V ) introduced in (13.21) still satisfies Proposition 13.31 namely 

di v B(u,V) = — Adiv V in D'(Q). (7.10) 

The proof of the above equality cannot be repeated as done in Section [5] because the minimizer u of 
the functional J P (Q) satisfies the regularity condition u £ IL 2 ’ 2 (n) only for p £ (1,3); nevertheless, 
for p > 3, one has that u £ W 2,g (fl) for any q < (p— l)/(p — 2) (see [15] Proposition 2.2]). Thus Vu 
is always in I / L 1 , 1 (H), and the distributional Hessian of u appearing in the definition B(u, V) is well- 
defined as a function in L 1 (H). Moreover, the function u satisfies the crucial regularity condition 
|V7t| p ~ 2 Vu £ W l ’ 2 (fX) (see [Tl] Corollary 2.1]), which enables us to obtain the proof of (17.101) as 
follows. Recall that cr = V/(Vu) = |Vu| p_ 2 Vu; then, since V/(-) is a locally Lipschitz function 
and Vu £ L°°(H) D VL 1 ’ 1 (H), we can apply the chain rule and we get Da = V 2 /(Vu)V 2 u a.e. in Q. 
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Thus the held B(u, V ) may be rewritten in terms of a as B(u , V ) = DaV — (DV — di vVI)a and 
in particular it belongs to L 2 (fl;M n ). It is then straightforward to obtain (17.101) . Namely, for any 
test function ip £ C£°(fi), it holds 

(div(DaV), ip) = -(DaV, Slip) = -{d k aiV k ,dnp) = (aid k V k ,dnp) + (aiV k ,d 2 k ip) 

= (— div(cf div V) + div(R div a) + di v(DVa), ip) 

= (div[(ZDI7 — div VI)a\, ip) — (A div V, <p ), 

where we have used the equality — diver = A (holding by (12.51) applied with g(t ) = —At). 

Now, we are in a position to rewrite m 0 0 similarly as in Theorem 13.41 (ii). Indeed, one can check 
that the proof of Theorem 13.41 (ii) given in Section O continues to work, thanks in particular to the 
identity (17.101) and to the W 2 ' q regularity of u for q < (p — l)/(p — 2)). Thus, defining Cd{u , V) as 
in (I3.10P (with f(z) = \z\ p /p), letting Q(u, •) be given by (17.71) . and setting By = (V, Vm), we obtain 


liminf^o r e(V) > rrioo = / Cd(u,V) d'H n inf Q(u,w — By) 

Jdn weCg°(Q) 


= / CD(u,V)dH n - inf Q(u,w — By ); 


(7.11) 


i an 


weHffin) 


here o(fi) denotes the space of functions in Hp’ 2 (fl) vanishing at dSl, and the last equality follows 
from the continuity of Q(u, •) on Wp' 2 . 

Step 2 (upper bound): Due to the growth of order p’ € (1,2] of the Fenchel conjugate f*(z*) = 
\z*\ p Ip', in order bound from above r e (V), we need to set up an approximation argument. Roughly 
speaking, we cover almost all D with a suitable family of subdomains of D in which the Hessian of /* 
is bounded from above, and then we pass to the limit. More precisely, as Vu is a continuous function 
which is strictly positive outside S, we can construct increasing sequence of open sets Qh t (^ \ S) 
such that, for every h £ N, |Vu| > 1/h in and the boundary := \ dSl is smooth and has 

a positive distance from dSl. For every fixed h, we follow the same procedure adopted for the proof 
of the inequality (14.81) in Lemma 14.21 in the variant when g is linear described in H7.ll We choose 
the fields r] k appearing in (17.31) as done in (17.41) , with the additional condition that p k = 0 in Q \ Q k 
(note that this is possible thanks to the assumption made on the support of the deformation field 
V). We infer that, for every fixed h, 


lirrisup r E (V) < m k '■= inf E*(ri,V ), 

£ 77£X(f2;M n ) , div rj= — X div V in f2 

77=0 in 


(7.12) 


where E*(rj,V ) is defined according to (17.11) (with f*(z*) = \z*\ p '/p'). 

By applying Lemma 18.51 and exploiting the regularity of u, we may rewrite rrih in primal form as 


m h = C D (u,V)d'H n inf Q h (u,w - By), (7-13) 

Jan w&v h 

where 14 denotes the subspace of H 1 (D/ l ) of functions having zero trace on dQ, and 

Qh(u, v) := 2 / (P(u)Vv,S7v) dx . 

J nh 
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Now, in order to bound from above the upper limit of rrih as h —> +oo, we are going to bound 
from below the infimum appearing at the r.h.s. of (17.131) . We start by noticing that such inhmum 
can be equivalently taken in the space Wj q (Q) of functions in Wp' 2 (yt) vanishing at dfl. Let 

{wh} C Wp q (11) be a sequence of minimizers for Qh(u,- — 6y). For every fixed k £ N, such 
sequence is bounded in so that it admits a subsequence weakly converging in to 

some element w^ k \ By a diagonalization argument, we can choose the same subsequence for every k, 
and hence the restriction of w^ to f li for l < k agrees with w®. Hence we obtain a (not relabeled) 
subsequence of Wh and an element w £ such that Wh converges weakly to w in iL 1 (H/ c ) for 

every k. Then, by Fatou’s Lemma and monotone convergence, we get 

lim inf Qh(u, Wh - 9y) > Q(u, w — 6y) > inf Q(u, w — 6y) . (7-14) 

h wffln) 

By combining (17.121) . (17.131) . and (17.141) . we obtain 

limsupr £ (H) < m*^ := f Cd(u,V) d'H n ^ 1 - inf Q(u,w — 9y). (7.15) 

e Jan weWp’g(n) 

Step 3 (conclusion): We finally need to show that the lower and upper bounds m oo and m ^ in (17.111) 
and (17.151) agree, and that they are equal to the expression at the r.h.s. of (17.91) . We firstly observe 
that, by arguing as in Remark 13.61 the integral over dQ of Cd(u,V) appearing in the expression 
(17.111) of moo can be rewritten as done in (17.91) . Then it only remains to prove that no Lavrenteev 
phenomenon occurs for the infimum problems appearing in (17.111) and (17.151) . In other words, after 
a translation, we are reduced to show that 

inf {Q(H, •) : w = —8y on dfl } = inf {Q(7q •) : w = —6y on dQ } . (7.16) 

w&wl' 2 (p.) w£h)’ 2 ( a) 

The validity of (I7.16|) is an immediate consequence of our assumption (17.81) and of the continuity of 
Q(u, •) in Wl' 2 {p). 

□ 


7.3 Perspectives 

A natural perspective is the possible extension of our results to the case of convex non-smooth 
integrands / and g. The main difficulty is that in this case one can no longer exploit the second order 
differentiability property for integral functionals stated in Lemma 18.11 However we expect that in 
some cases the second order shape derivative still exists, and admits a representation formula similar 
to (13.61) . with the quadratic form Q(u, •) (formerly associated with the second order differentials of 
/ and g ) replaced by a suitable convex, positively 2-homogeneous, non-quadratic function. This 
kind of generalization would be useful in applications: for instance, the shape functional obtained 
by taking the function f(z ) equal to ^ if |^| > 1 and |z| if \z\ < 1, and the function g(v) 

equal to —Xv (for some A £ M), is related to the optimization of thin rods in torsion regime (see 
Eli), and studying its second order shape derivative might be helpful in order to investigate the 
occurrence of homogenization regions in an optimal design. (For the computation of the first order 
shape derivative in this case, see Example 3.9 (ii) in |8j.) 

Another perspective is trying to understand whether the results obtained in (17.21 may entail some 
useful information on the second order shape derivative of the shape functional Joo(H) obtained as 
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the limit as p —> +oo of the p-torsional rigidity functionals J p (fi) defined in (13.161) . Actually, up to 
constant multiple, it is well known that here holds Joo(fl) = Jq dgn(x) dx and lim p u p (x) = dgn(x), 
where dgQ is the distance function from dLl, and u p is the unique solution to J p (Ll) (see @J[7], (57)). 


8 Appendix 

Lemma 8.1. Let <f> : R n —» R be a strongly convex function of class C 2 satisfying growth conditions 
of order 2 from above and below, and let 1^ be the integral functional defined on L 2 (P; R n ) by 

■= [ <j>(z(x))dx. (8.1) 

Jn 

Let zq be a fixed vector field in L°° (LI] W 1 ) and, for h £ L 2 (P; R n ), setq^(h) := ^ f n (V 2 <j)(zo)h, h) dx, 
and 

A ^ h( z o + eh) - I^(zp) - e(V^), h) L i ^ ^ 

Then: 

(i) Ig is of class C 1 on the Hilbert space L 2 (Ll] R n ), and XL^(z) = V(f(z); 

(ii) A e ^ Mosco converges in L 2 (P;R n ) to qg as £ —> 0; 

(in) L'f, is second order differentiable at zo with respect to L°° variations; namely, for every h £ 
L°°(P;R n ) ; it holds A £ ^(h) —> q^h) as e ^ 0. 

Proof. For statements (i) and (ii), see [351 Proposition 4.1 and Corollary 3.2], Let us now prove 

(iii). For a fixed h £ L°°(P;R n ), the Mosco convergence property stated in (ii) implies that 

q<b(ti) < lim inf A e Ah). 

£->0 

In order to prove that limsup e A £ ,(j>(h) < q^(h), we exploit the second property given by the Mosco 

convergence: there exists a sequence h £ £ L 2 (P; R n ) such that h £ —>• h and q$(h £ ) = lim £ _>.o A e ^(h e ). 
We remark that the function A £j ^ is convex, therefore, A £ ^(h) — A £t< j,(h £ ) < (VA £t <f,(h),h — h £ ) L 2 . 
If we prove that, for e > 0 small enough, 


\\VA £ ^(h)\\ L2 <C, 


(8.3) 


we are done, since we infer 

lim sup A £ ^(h) - q^h) = lim sup (A £ ^(h) - A £ ^(h e )) < lim C\\h £ - h\\ L 2 = 0. 

We conclude by proving the claim (18.31) : we remark that VA £ (h) = [\7(j)(zo + eh) — X 4 >(zq)]/ e, and 
hence ([8.3D follows from the local boundedness of the Hessian matrix X 2 (j) and recalling that by 
assumption there exists R > 0 such that \\zq + A R for t small enough. □ 


Lemma 8.2. Let X be a normed vector space and let X* be its topological dual. We denote by 
(■, ■) the duality product. Let h : X — > R U {+ 00 } be a proper function. Then, for every a £ X and 
b £ X*, there holds 

(h(- - a))* (z*) = h*(z*) + (a, z*) , (h(-) + ( b , ■»* (z*) = h*(z * - b). 
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Proof. The statement follows straightforward by using the definition of Fenchel conjugate. □ 

Lemma 8.3. Let A be a positive definite n x n matrix and let Qa denote the associated quadratic 
form, defined as Qa(x) := \(Ax,x). Then (Qa)* = Qa- 1 - 

Proof. See [38i Chapter III], □ 

Lemma 8.4. Let X be a Banach space and let h : X -A- M be a strongly convex function of class 
C 2 . Then dom(h*) has nonempty interior, h* is C 2 on the interior of dom(h*), and x* = Xh(x ) 
with x € X implies x = Vh*(x*) and X 2 h*(x*) = ( V 2 h(x )) 1 . 

Proof. See m Proposition 10 in Section 2 of Chapter II]. □ 

Lemma 8.5. Let Y,Z be Banach spaces. Let A : Y —>■ Z be a linear operator with dense domain 
D(A). Let 4> : Y -A- MU{+oo} be convex, and T : Z —> MU{+oo} be convex lower semicontinuous. 
Assume there exists no € D(A ) such that <h(tto) < +oo and T is continuous at Auq. Let Z* denote 
the dual space of Z, A* the adjoint operator of A, and 4>* ; T* the Fenchel conjugates of&, T. Then 

— inf |4 r (Au) + 4>(u)| = inf j\H*(<r) + <!>*(— A* <r)j , (8-4) 

and the infimum at the right hand side is achieved. 

Furthermore, u and a are optimal for the l.h.s. and the r.h.s. of m respectively, if and only if 
there holds a G ddt(Au) and — A*a € (9<h(lt). 

Proof. See [BJ Proposition 14]. □ 
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